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We show that Popov-Fedotov trick of mapping spin-1/2 lattice systems on two-component 
fermions with imaginary chemical potential readily generalizes to bosons with a fixed (but not lim- 
ited) maximal site occupation number, as well as to fermionic Hamiltonians with various constraints 
on the site Fock states. In a general case, the mapping — fermionization — is on multi-component 
fermions with many-body non-Hermitian interactions. Additionally, the fermionization approach 
allows one to convert large many-body couplings into single-particle energies, rendering the dia- 
grammatic series free of large expansion parameters; the latter is essential for the efficiency and 
convergence of the diagrammatic Monte Carlo method. 
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It has been demonstrated recently that Feynman's di- 
agrammatic technique can be used for controllable first- 
principles numeric studies of lattice fermions (the Hub- 
bard model at U/t < 4) [1]. From the diagrammatic 
point of view, lattice fermions feature a unique combina- 
tion of three properties: (i) Wick's theorem for averages 
over the non-interacting system which allows to expand 
results for the interacting system into Feynman's dia- 
grammatic series, (ii) restricted Hilbert space prevent- 
ing the system from the Dyson's collapse argument [2] 
thereby allowing for the possibility that the diagram- 
matic series has a finite convergence radius, (iii) sign- 
alternation of diagrams within the same expansion order 
leading to the so-called 'sign-blessing' phenomenon when 
high-order diagrams cancel each other with factorial ac- 
curacy. Spin systems (as well as bosonic ones with re- 
stricted filling factor) share the property (ii), but, speak- 
ing generally, not the properties (i) and (iii). Bosons 
with unrestricted filling suffer from the Dyson's collapse. 
An idea of how the problem of Dyson's collapse could 
be circumvented for classical fields and bosons has been 
proposed recently [3] ; the price for the finite convergence 
radius of the regularized theory is complexity of the inter- 
action vertices which increases with decreasing the (con- 
trollable) systematic error. An outstanding problem of 
numerics based on Feynman diagrams, even with for- 
mally convergent series, is the regime of large coupling 
constants, where the bare perturbative series fails to con- 
verge at small enough orders of expansion amenable to 
treatment by the Diagrammatic Monte Carlo (DiagMC) 
method [1]. 

In this report, we show how the Popov-Fedotov trick 
[4] of representing spin-1/2 lattice systems with two- 
component fermions having imaginary chemical poten- 
tial and its extensions to other spin values and projected 
Hamiltonians [5], can be generalized to arbitrary spin, 
boson and fermion models with various constraints on 
the site Fock states. We would like to stress here the dif- 
ference in the final formulation of the problem between 



the Popov-Fedotov type theories and the slave fermion 
approach, see e.g. [6]. While constraints on the Hilbert 
space in the latter method are implemented with the use 
of additional local gauge fields, in Popov-Fedotov's (and 
our) formulation the same goal is achieved by construct- 
ing an equivalent Hamiltonian based solely on conven- 
tional fermionic particles. In the general theory presented 
below all constraints on the site Fock states are 'relo- 
cated' from the Hilbert space into the system's Hamilto- 
nian by introducing non-Hermitian terms with imaginary 
couplings. The power of this approach is not limited to 
providing universal tools for constructing the Feynman's 
diagrammatic technique with presumably finite conver- 
gence radius for arbitrary spin/boson systems. Yet an- 
other crucial feature of the method is the possibility of 
absorbing large on-site coupling constants of the origi- 
nal model into single-particle energies of additional ficti- 
tious fermions. The remaining interactions between the 
fermions are either of the order of unity (in terms of the 
hopping amplitudes), or scale linearly with temperature. 
This feature can prove useful for treating purely fermionic 
systems — e.g., the Hubbard model with large on-site in- 
teraction. 

Fermionization. Let us start with bosons and spins. 
Here the idea is to represent bosonic/spin degrees of free- 
dom by internal states of fictitious fermions living on 
the sites of the corresponding lattice. To ensure that 
the physical Hamiltonian matrix is not modified by the 
fermionic nature of fictitious fermions the new Hamilto- 
nian does not have terms involving hopping of particles. 
All matrix elements of the original Hamiltonian are then 
encoded in interaction processes which change internal 
degrees of freedom of dispersionless fermions immediately 
leading to the dynamical equivalence of the original and 
'fermionized' Hamiltonians by construction. The prob- 
lem arises with the statistical equivalence. The Hilbert 
space of the fermionized Hamiltonian, H, contains non- 
physical Fock states of the sites (e.g. with occupation 
numbers other than unity) and thus its partition func- 
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tion differs from that of the original Hamiltonian. The 
trick is to introduce additional terms, H*, which, on one 
hand, produce zero when acting on physical Fock states 
of the sites (and thus are dynamically irrelevant for the 
physical subspace), while, on the other hand, make the 
partition function of H' = H + identical to that of 
the original Hamiltonian. The crucial observation is that 
this goal is readily achieved with involving only on- 
site interactions, H„ = ^ . (here j is the site index), 
due to the following structure of the trace of H' 

where Z is the physical partition function (sum over all 
physical states of all the sites), £aA labels different con- 
figurations of J\f non-physical sites, and Z^ M is the parti- 
tion function of a physical system living on a lattice from 
which all the non-physical sites are removed. Finally, Sj 
is a single-site trace over all non-physical states of the 
site j. The generalization of Popov- Fedotov trick is to 
chose in such a way that Sj — 0. We will see that 
there exists a continuum of straightforward options for 
doing that individually for each site j. 

Obviously, the idea of relocating constraints on the site 
Fock states from the Hilbcrt space into the Hamiltonian 
can be used in a more general context. A characteris- 
tic example is provided by the t-J model dealing with 
two fermionic components subject to the constraint that 
there are no doubly occupied sites. Since there is only 
one non-physical Fock state per site, one cannot achieve 
the situation of Eq. (3) with Sj = 0. The trick is to 
introduce an auxiliary on-site fermionic mode that does 
not interact with the rest of the system and has zero 
Hamiltonian in the physical sector, i.e. it leads to trivial 
replacement Z —> 2 V Z, where V is the total number of 
system sites. For non-physical sites, however, coupling to 
the auxiliary mode is chosen in such a way that Sj = 0. 
We find it reasonable to use the term fermionization in 
this case as well, since the formalism of constructing the 
Hamiltonian H' here is often the same. Since introduc- 
ing an auxiliary mode is not subject to any restrictions, 
the trick is universal. 

Formalism. Consider first the case involving the above- 
mentioned auxiliary fcrmion mode. One solution is based 
on Eq. (3), meaning that we need the Hamiltonian H 
to be written in fermionic operators in such a way that 
the physical sector is an invariant subspace of H. This 
is readily achieved by multiplying each elementary term 
in H by projectors onto the physical states of all the 
sites this term deals with. Note that for some terms 
this procedure is redundant; an important example being 
terms diagonal in the site Fock representation. Then we 
take 

H { J } = inTP^hj , (2) 



where P* is the projector onto non-physical Fock states 
of the site j, and hj is the number operator for the auxil- 
iary fermionic mode. By construction, H and com- 
mute, which brings us to Eq. (3). In view of (2), we then 
have Sj oc Yla=o l exp(— wrri) = 0. 

Less restrictively, one can demand that H is only pre- 
serving the notion of a non-physical site, i.e. it does not 
change the £v configuration but is allowed to change non- 
physical states of the sites and couple them to physical 
ones. Now, in the expression 

Tre-«'/T = 2 v z + £ 2 V ~ N Z^ E *i > ( 3 ) 
A/"=i Or je?Ar 

is the partition of H, including non-physical sites, 
and excluding the auxiliary-fermion sector — the latter is 
evaluated explicitly due to factorization — for a given set 
£/V of non-physical sites; Sj = J2n=o l exp(— i7m) = 0. 

Now we turn to a more specific case of spins/bosons 
(which allows to avoid the auxiliary fermion trick) and 
consider the fermionized Hamiltonian H which includes 
both on-site, Hq, and inter-site, Hi, parts: 

H = H + Hi , (4) 

The inter-site Hamiltonian has the following structure 

^ = E E *Zp„*Q%+fiQ%- (6) 

i^j a,/3,7,(5 

Here i and j are the site labels, Greek letters stand for 
the physically relevant fermionic Fock states of a site, 
and ^ s is encoding all non-local matrix elements of 

the original physical system. The on-site operator Q^\g 
changes the state of the site i from a to /3, yielding iden- 
tical zero otherwise (most notably, for all the physically 
irrelevant states of the site) . With the diagrammatic ex- 
pansion in mind, we need this operator to be a polyno- 
mial in terms of the fermionic creation and annihilation 
operators. That can be guaranteed by writing it in the 
form 

q% = . ( 6 ) 

where A^_^g is the polynomial of the fermionic creation 
and annihilation operators, such that when acting on the 
state a it produces the state (3, and is a projector 
on the state a. 

For example, if we want to 'fermionize' a system of 
single-component lattice bosons (with the constraint that 
the maximal site occupancy is mo), the physical state of 
a site with m < mo bosons is [6] (here we suppress the 
site label for clarity) 

| a = m) = | n m = 1; rii = 0, I ^ m ) . 

Here n; = 0, 1 with < I < rriQ is the site occupa- 
tion number for the l-th fermionic component. Since 
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physically relevant states have exactly one fermion per 
site, for single-particle hopping transitions we have (be- 
low I = m ± 1) 
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(7) 



where / t - m (fi m ) is the creation (annihilation) operator 
for the m-th fermionic component on the site i, and = 
fik fik- [There is no need to include n, m (l— riu) operators 
in the state projector since matrix elements of are 
automatically zero unless m m = 1 and nu =0.] 

The on-site Hamiltonian — for simplicity, we discuss 
here the case when Ho is diagonal in the site Fock basis — 
can be cast into different forms. By insisting that H is 
zero for non-physical sites we arrive at 



(8) 



where E a is a real number. As long as all E a are not 
large compared to the non-local matrix elements, this 
form appears to be sufficient. When some parameters 
E a become large (note that for bosons E a describe the 
on-site interaction), the form (8) becomes problematic 
for the diagrammatic technique because its non-linearity 
in fermionic occupation numbers entails the necessity to 
expand in powers of E a . It turns out that in such cases 
one can switch to Hq where projectors are replaced with 
occupation numbers of fermionic modes, if necessary, at 
a price of adding auxiliary degrees of freedom. We will 
start with Hq in the form (8), and then upgrade the the- 
ory to the general case. 

Consider the trace of e~ H l T in the site Fock represen- 
tation. Since all terms in H contain projectors on the 
physical states, it gives zero when acting on any non- 
physical state of the site. Likewise, by construction all 
matrix elements of are zero for each physical site 
meaning that [H, H*\ = 0. As a result, we arrive at the 
structure of Tre~ H l T expressed by Eq. (3). Indeed, for 
a given configuration (A/", £v) of the- non-physical sites, 
we can ignore these sites when calculating the trace of 
£ -h/t ^ anc j jg n0 re physical sites when calculating the 
trace of e~ H */ T '. This leads to the factorization of the 
partition function, Z(H',£jy) = Z(H,£tf) Z(H*,£m), 
since traces in Z(H,£jy) and Z(iJ*,£/v) = L] 
taken over mutually exclusive sets of sites. 

To have Sj = 0, one can choose if* in the form (since 
H+' acts only on non-physical states it does not need to 
be Hcrmitian) 
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(9) 



where a labels non-physical states of the site j, and w & 
are certain complex numbers. With Eq. (9) we have Sj — 
S = Wa and to satisfy the requirement S = it is 



necessary and sufficient to choose complex numbers w 6 
in such a way that their sum is zero, ^ a w & = 0. 



One simple possibility is to consider u> a 
tune phases 
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(10) 



Consider, e.g., the Popov-Fedotov case of spin- 1/2 sys- 
tem. The two physical states arc |1) = \n^ = l,nj, = 0) 
and |2) = = 0,n^ = 1). The two non-physical states 
are |3) = |n-|- = 1,714. = f) ana - 1^) = l n t = 0j n i = 0)' 
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■rt n tt and P i ] = (f - "jtX 1 - ^4-)- Wc 
see that the original Popov-Fedotov trick corresponds to 

03 = -04 = ±7r/2: 
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Obviously, there is a continuum of other choices. For 
example, 03 = it, 04 = yielding if* = —iTirnj^riji. 

Eliminating large on-site coupling. In fermionized 
spin/boson systems one can always avoid the problem 
of large coupling constants associated with large values 
of E a . Consider lattice bosons with mo > 1 where strong 
on-site interaction U translates into large values of E a . 
In this case, we chose Ho in the diagonal bilinear form of 
fermionic operators: 



m 



(if) 



Though Ho docs not nullify the non-physical states, it 
still commutes with the purely local H* and preserves 
factorization of the partition function in Eq. (3). Wc 
simply have to account for Ho energies when calculating 
iSj-factors. For H* given by Eq. (9) we then get 



-V&/T, 
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E m (a\n m \a) 



(12) 



(13) 



We see that achieving S = is easy if for any state 
I a) there exists at least one other state \a') with ex- 
actly the same value of rj & — states with the same energy 
can be made canceling each other by using phase factors 
Wa = e 1 ^ 3 only. [Needless to say that the requirement 
of having two non-physical site states with exactly the 
same energy can be always met by simply introducing 
the abovc-discusscd auxiliary fermionic mode] 

Less obvious is the fact that there always exists a 
Popov-Fedotov type solution for a generic bosonic/spin 
problem with the Hamiltonian (4)-(5), (11) and local 
H { J ] of the form [7] 
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Indeed, Hi 3 ' commutes with the physical Hamiltonian 
and is identically zero in the physical sector. The expres- 
sion for Sj then has the form 

S J =V ma (X 1 ), (15) 

where 'P mo (Xj) is an order-mo polynomial in Xj. By the 
central theorem of algebra there exist mo (including de- 
generacies) solutions to (15). The choice of Xj does not 
matter; one is free to choose one value for all sites or 
work with their arbitrary site distribution. 

Second fermionization. Fermionic Hubbard model. 
The fermionization technique for bosonic/spin modes can 
be also used in the context of purely fermionic sys- 
tems/modes to eliminate direct on-site interactions be- 
tween the fcrmions. The general principle is to exclude 
from the system's Hilbert space states with more than 
one fermion on a site by interpreting them as either hard- 
core composite bosons (for an even number of original 
fcrmions on a site), or composite fcrmions (for larger than 
unity odd number of original fermions on a site). The 
composite nature of the new bosonic/fcrmionic modes 
is not relevant for the final theory which treats bare 
and composite particles on equal footing. The hard-core 
bosonic modes are then fermionized; we call this proce- 
dure second fermionization. 

As an example, consider the fermionic Hubbard model 

H = - ^2 fjvfio + \. Un n n u ~ mt'M ~ H n i\\ > 

<ij>,a i 

where a ="f , -J-, and < ij > restricts the summation to the 
nearest-neighbor sites (hopping amplitude is set equal to 
unity) . We are interested in the case when the on-site in- 
teraction U and the chemical potentials fi a can be large. 
We interpret the doubly occupied sites, = = 1, as 
containing one hard-core boson (doublon) with the en- 
ergy = U — — fi±. Within this representation, 
the Hamiltonian is re-written in terms of the fermions 
subject to the constraint + < 1 and the doublons 



subject to the constraint mo = 1. The constraints on 
fcrmions and doublons — upon fermionization of the lat- 
ter using fictitious two-component fcrmions — are imple- 
mented by using tools of projectied operators to deal with 
hopping matrix elements, Eqs. (5), (6). The elimination 
of non-physical sites can then be achieved by the auxil- 
iary fermionic mode. The resulting model thus will in- 
volve five, instead of original two, fermionic modes. How- 
ever, in the new Hamiltonian all multi-particle coupling 
constants are equal to unity (the imaginary coupling con- 
stant in iJ* scales as T). 
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